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Scale-free networks with topology-dependent interactions are studied. It is shown that the uni- 
versality classes of critical behavior, which conventionally depend only on topology, can also be 
explored by tuning the interactions. A mapping, 7' = (7 — (J,)/(l — fi), describes how a shift of the 
standard exponent 7 of the degree distribution P(q) can absorb the effect of degree-dependent pair 
interactions Jij oc (ftg'j) - ''. Replica technique, cavity method and Monte Carlo simulation support 
the physical picture suggested by Landau theory for the critical exponents and by the Bethe-Peierls 
approximation for the critical temperature. The equivalence of topology and interaction holds for 
equilibrium and non-equilibrium systems, and is illustrated with interdisciplinary applications. 
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In this Letter we pose and answer the following funda- 
mental questions. What are the relevant variables that 
determine the universality classes of critical behavior in 
networks? Is the conventional classification based on the 
exponent 7 of the distribution function of network con- 
nections P(q) sufficient? What is the effect on the critical 
behavior of interactions that depend on the connectivity? 
Do the interactions give rise to new relevant variables or 
can these be "transformed away" topologically? We fo- 
cus on scaling arguments and analytic approaches 

Nowadays, scale-free networks enjoy a lot of attention 
due to their ubiquitous occurrence and the suitability 
of modern computational techniques for understanding 
their properties and predicting their behavior 0. These 
networks typically have simple topology. The usual no- 
tions of spatial coordinates and dimensionality are not 
present; e.g., while the "volume" scales as the number of 
nodes N, the largest distance or "diameter" grows not 
faster than logiV, suggesting infinite "dimension". This 
is why mean-field or Landau theories are successful. 

Understanding critical phenomena in scale-free net- 
works is of fundamental and practical relevance. There is 
a surprising diversity of "mean-field" universality classes, 
which are topology-dependent 3]. However, for most real 
networks and for the simplest interactions between nodes, 
the critical point is inaccessible (e.g., T c = 00). We pro- 
pose a way around this by tuning the interactions. 

Without loss of generality, consider first the context 
of opinion formation in "sociophysics" 4J. The network 
consists of people. Opinions are represented by spin ori- 
entations, and communication by pair interactions. The 
statistical physics of the formation of a common opinion 
is akin to that of spontaneous symmetry breaking be- 
low a critical temperature T c in a Hamiltonian model for 
equilibrium cooperative phenomena p|. 

Putting Ising (or Potts) spins s, on the nodes i = 



1,...,N and ferromagnetic interactions J on the edges, 
leads to interesting critical behavior, different from its 
counterpart in lattice spin models For the standard 

scale-free distribution P(q) oc q~ 1 of the number of con- 
nections or "degree" of a node, the decay exponent 7 is 
the key parameter distinguishing the universality classes. 
Standard mean-field critical behavior (a = for the spe- 
cific heat, (3 = 1/2 for the order parameter, j sus = 1 for 
the susceptibility) is predicted for 7 > 5. Logarithmic 
corrections appear at 7 = 5. For 3 < 7 < 5 the critical 
exponents a and (3 are non-universal; they depend on 7. 
The susceptibility exponent 7 stis = 1 is superuniversal. 
Conventional finite-size scaling (as a function of network 
size N) predicts the familiar constraint a + 2f3 + ^ sus = 2. 
For 7 < 3, however, the framework of finite-temperature 
critical phenomena breaks down as T c moves to infinity. 

Which values of 7 apply to real scale-free networks? 
Most studied ones, including WWW, Internet, collabo- 
ration, citation, cellular, ecological or linguistic networks 
have 2 < 7 < 3. The Barabasi-Albert (BA) network, 
grown by preferential attachment, has 7 = 3. 

The initial motivation for introducing degree- 
dependent interactions Jij = J(qi,qj) was to prevent 
T c from diverging in the Ising model on the BA net- 
work [j| by compensating high connectivity with weak 
interaction through Jy oc 1/ \Miq~j |7|. Further work has 
led us to observe that this interaction is equivalent to 
a g-independent J, provided 7 is shifted from 3 to 5. 
A Landau-theoretic argument and a mean-field ansatz 
generalize this equivalence (see further) and provide the 
insight that interactions that depend on connectivity 
effectively modify the topological distribution P(q). 

We introduce a family of interactions parameterized by 
an exponent /x, 
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where (.) denotes the average over the degree distribution 
P(q). Generalizing the Landau theory of Goltsev et al. 
we obtain the following form, in zero external field, for 
the "free energy" $ as a function of the order parameter 
x = J^i Er( 8 i)/N, with Et(-) the thermal average, 

$(x) = (<P(q-»x,q 1 ->*x)), (2) 

where the first argument of the constrained free energy 
<j> is the rescaled order parameter, and the second one is 
the rescaled effective field acting on a node with q con- 
nections. Both rescalings, by a factor come from the 
qi- and qj-dependence of Jy. The crucial assumption Q 
is that <p(y, z) can be written as a power series in y and 
z. The singularity structure of $>(x), leading to fasci- 
nating deviations from standard mean-field behavior, is 
then induced by the fact that the moments (q n ) diverge 
for n > n c =7 — 1. This in turn leads to a divergence 
of the coefficient f n of x n in the free energy. Thus, in 
this system non-classical critical behavior is possible, not 
due to spatial fluctuations of the order parameter (since 
space is effectively infinite-dimensional), but due to the 
scale-free character of the degree distribution. 

Our main result is that a network with exponents (7, fi) 
can be mapped onto one with (7',// = 0), in the sense 
that both are in the same universality class of critical 
behavior. The latter has constant couplings, independent 
of q. We find the exponent relation 

Y = (7-M)/(l-A0, (3) 

which can be proven as follows. The critical value n = n c 
for diverging coefficients /„ in the free energy must be the 
same in order for the two models to be equivalent. Since 
the leading moment (q n ) in /„ in the model with fi = 
gets replaced by (c^ 1- ^)™) in the model with /1 ^ 0, we 
obtain n c = Y — 1 = (7 — 1)/(1 — fi). 

A more general "mean-field" proof, which does not 
rely on a free energy, runs as follows. Within a mean- 
field approach the only way in which the degree q enters 
in physical properties is through the quenched average 
interaction, over all networks, between any two nodes 
with fixed degrees qi and qj . This average is JijPij , with 
Pij = QiQj/((q)N) the probability that i and j are con- 
nected. Note that even for constant J (case fj, = 0) the 
quenched average interaction is (/-dependent, while for 
/i = 1 it is not. For /i / 0, the qi are transformed 
to q[ = q^ , using (1). In order to retain the same 
physics, averages over the degree distribution must be 
invariant. This requires a distribution transformation, 

P(q)=P'(q'(q))dq'(q)/dq, (4) 

from which (3) follows for scale-free (power-law) P(q). 

Consequently, for a scale-free network the range \x £ 
[2 — 7, 1] allows one to explore the whole range of uni- 
versality classes uncovered in previous works. It is no 
longer necessary to vary the network topology. It suf- 
fices, for a fixed 7, to tune the form of the interaction. 



An important consequence is that for real scale-free net- 
works, with typically 7 < 3, one is no longer set back 
by an infinite critical temperature when putting on in- 
teractions. For the BA network the limit /i — > 1 maps 
onto 7' — » 00 and corresponds to the cross-over to a thin- 
tailed degree distribution of, e.g., Poisson type. In the 
opposite limit, fj, — » — 1, for which 7' — > 2, the threshold 
is reached beyond which the first moment (q) diverges in 
the equivalent network with constant couplings. 

For locating the critical temperature a Bethe-Peierls 
(BP) approximation, which emphasizes the local Cayley 
tree-like structure of the network, normally gives a very 
reasonable first approximation [(|. The entropic contri- 
bution to the free energy is truncated to single-spin and 
pair terms @. Besides this standard assumption, we im- 
pose a new scaling relation on the local order parameter, 
Et{s{), in a given network. This scaling is in harmony 
with the definition of the effective field acting on Si in 
the Landau theory and is corroborated by MCS. It reads 

E T ( Sl )^q]-»x/{q^). (5) 

Using this to eliminate Sj in terms of qi and x one obtains 
the self-consistent BP equation near T c , 

(g)^P(g)( g -l)^ = 

^P(gO^E P fe)? 2 (l + tanh^M!!_) \ (6) 

qi q-2 ^ B 1 2 ' 

For the conventional case \i = 0, we can extract from 
this the large-(g) approximation kBT c / J « (q 2 )/(q) — 1, 
in agreement with exact results @. 

We now turn to illustrations of this framework and 
start with 7 = 3 and /i = 1/2 ("Special Attention Net- 
work" 0). We discuss the critical point, specific heat, 
order parameter and susceptibility. 

The critical temperature versus average degree (q) is 
shown in Fig.l, for (q) = 2, 4, and 6. It has been derived 




FIG. 1: T c versus (q) for a scale-free network with 7 = 3 and 
H — 1/2. BP estimates (filled squares) almost coincide with 
exact results (replica technique, open stars). The solid line is 
the lattice BP approximation for Q = (q) and the dashed line 
is the MF approximation. 
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by exact solution of the model using the replica tech- 
nique. Since all couplings are "ferromagnetic" the replica 
symmetry is not broken. For comparison, the almost co- 
incident results from the BP approximation (6) are also 
indicated. A good rule of thumb is ksT c /J ps (q) — 1. 
Incidentally, this is also the result of the large- (q) expan- 
sion which can be derived from (6). Interestingly, the 
results are very close to the conventional BP approxima- 
tion J/k B T c = 0.51n[(Q - 2)/Q] (thin solid line), for a 
regular lattice with coordination number Q, except for 
Q = 2. For completeness the mean- field conjecture Q 
kBT c /J = (q) is also shown. 

These calculations assume no correlations exist be- 
tween the edges of the network. The rule of preferen- 
tial attachment violates this assumption, so the BA net- 
work must be considered separately. For the special case 
(q) = 2, a BA network differs strongly from a correlation- 
free network. The former consists of a single tree-like 
structure without loops, whence T c = 0. In contrast, an 
uncorrelated network with (q) — 2 can consist of clusters, 
which can have loops, leading to a finite T c , as replica 
technique and BP approximation predict. 

Having located the critical temperature, and eluci- 
dated its non-universality through its dependence on (q) 
and on the network correlations, we now turn to more 
universal properties. The specific heat singularity for the 
Ising model on scale-free networks has a very interesting 
and subtle form 6]. For fi = 0, it varies from a classic 
mean-field jump, predicted for 7 > 5, to a continuous be- 
havior but with diverging slope, for 4 < 7 < 5, and with 
continuous slope for 3 < 7 < 4. For // = 1/2, the equiv- 
alence relation, Eq.(3), predicts, with r = (T — T c )/T c , 



C s 



cx (lnr -1 )" 



(7) 



The shape of this specific heat near T c is astonishingly 
different from that of the conventional mean-field jump. 
Although the critical exponent a consistent with (7) is 
zero, as for standard mean field, the inverse logarithm 
ensures a vanishing jump, with diverging slope. 

Evidence gathered from MCS, and, more convincingly, 
the cavity method, shows that i) the specific heat C(T) 
reaches a maximum well below T c and ii) the jump sin- 
gularity that is apparent for small network size N closes 
slowly when TV is increased. These findings are consistent 
with (3), and suggest that for 7 = 3 and fi — 1/2, the 
network indeed maps onto one with 7 = 5 and // = 0. 

Figures 2 and 3 illustrate these results for the spe- 
cific heat, for a BA network with fi = 1/2. Shown are 
MCS (Fig.2) for (q) = 10 and N = 5600, and the cav- 
ity method (Fig. 3) for (q) — 4, applied to the sequence 
N = 100, ...,10 6 . The latter technique clearly suggests 
the slow closing of the specific heat jump, for N — > 00. 

We have also computed the order parameter ("mag- 
netization") and susceptibility singularities near T c and 
find, for \i = 1/2, that the results are consistent with 
j3 = 1/2 and j svs = 1. Again, these values agree with 
the predictions from the mapping (3). It should be re- 
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FIG. 2: Specific heat from MCS of BA networks (7 = 3), 
with fi = 1/2 and (q) = 10. The critical temperature in the 
large- N limit is k B T c /J « 8.96 (Eq.6). For N = 5600, MCS 
of the susceptibility maximum indicates fcsT c (A r )/J w 8.65, 
well above the temperature of the specific heat maximum. 
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FIG. 3: Specific heat from the cavity method for BA net- 
works, with fi = 1/2 and (q) = 4. Shown are data for N = 10 2 
(squares), N = 10 4 (triangles), and TV = 10 6 (stars). Lines 
are guides to the eye. The finite-size critical temperature (ap- 
parent from the jump in C) rapidly converges for large N, to 
ksTc/J ~ 2.92 (Eq.6). Clearly, the maximum in C precedes 
T c and the jump in C tends to close for large N. 



marked that for 7 = 5 and /i = 0a logarithmic correction 
factor is predicted for the order parameter, which, how- 
ever, cannot easily be detected on top of the square-root 
singularity (3 = 1/2. Further, our result for j sus is not 
discriminative, since jsus = 1 is superuniversal, valid for 
all 7 > 3 in the model with // = 0. 

Further evidence for (3) is obtained for the interesting 
case 7 = 3 and /i = 1/3. A linear specific heat (cv = 
— 1) and a linear order parameter singularity (/? = 1) 
clearly emerge for large N in cavity method computations 
(Fig. 4). MCS of specific heat (Fig. 5), order parameter 
and susceptibility supports this conclusion. The results 
agree with what is expected for 7' = 4 and /a' = 0. 

Additional evidence for (3) comes from a study of net- 
works with 7 = 3 and ji = 1. Using MCS and cavity 
method the standard mean-field jump of the specific heat 
is retrieved. This is consistent with a shift of 7 to a value 
greater than 5 in the reference system with /i = 0. 
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FIG. 4: Specific heat from the cavity method for BA net- 
works, with fi — 1/3 and (q) = 4, for N — 10 2 (circles), 
N = 10 3 (squares), and N = 10 5 (stars). A linear In- 
dependence develops for large N and meets the high-T back- 
ground with a jump in slope at T c . The finite-size critical 
temperature converges to fesT c /J ~ 4.21 (Eq.6). 



We now leave equilibrium statistical mechanics and fo- 
cus on dynamical systems. In the contact process for 
disease spreading each node can be either ill or healthy. 
An ill node can cure at rate 1, and a healthy node be- 
comes infected at a rate which is A times the number of 
ill neighbors. The model has a phase transition between 
an absorbing healthy state and an active state with a 
non-zero density of ill nodes, at some A c . On a scale- 
free network a finite A c is found for 7 > 3. The critical 
exponents depend on 7 in the range 3 < 7 < 4 and 
assume standard mean- field values for 7 > 4 [9(- Gener- 



alizing the process to a degree-dependent infection rate, 
Ay = A(g) 2/i (gj(7j)~' 1 , we arrive again at (3). The exp o 
nent mapping appears to be quite generally valid [l 
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FIG. 5: Specific heat from MCS of BA networks with N = 
5600, for /1 = 1/3 and (q) = 10. The critical temperature in 
the large- N limit is at k B T c /J « 11.33 (Eq.6). 

In conclusion, static and dynamic order-disorder tran- 
sitions on scale-free networks display singularities that 
depend on the network topology and on the form of the 
interactions. Connectivity-dependent interactions can be 
used as a probe of topology-dependent cooperative be- 
havior. The exponent mapping (3) prescribes how to 
trade interactions for topology. 
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